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Accompli  shments 

The  main  accomplishments  of  the  research  carried  out  under  this  contract 
are  discussed  in  this  section.  The  first  year  was  devoted  to  computer  coding 
the  rotational  Hamiltonians  for  the  various  vibrational  states  in  methane. 

The  second  and  third  years  were  devoted  to  a  new  formulation  of  the  model 
Hamiltonian  based  on  a  jellium  approximation  for  the  molecule.  Finally, 
in  the  third  year,  consideration  was  given  to  a  calculation  of  the  vibrational 
modes  of  simple  molecules  based  on  the  electronic  structure  of  these  molecules. 
This  latter  work  is  continuing  under  AFOSR  grant  #77-3130. 

The  calculation  of  the  rotational  levels  of  the  high  symmetry  molecules 
was  carried  out  using  standard  techniques  developed  for  simpler  low  symmetry 
molecules.  The  computer  code  was  written  with  the  help  of  Samuel  Safran 
who  finished  a  Ph.D.  degree  at  M.I.T.  and  then  obtained  a  position  at  Bell 
Laboratories.  It  was  shown’ that  the  Wang  representation  for  the  angular 
momentum  matrices  gives  rise  to  a  factorization  of  the  rotational  Hamiltonian 
for  all  symmetrices  of  vibrational  states.  This  general  result  is  applied 
to  the  methane  molecule  below. 

The  molecular  vibration  rotation  energy  levels  for  the  methane  molecule 
are  formulated  making  full  use  of  the  high  symmetry  of  this  molecule.  The 
vibration-rotation  Hamiltonian  in  Td  symmetry  obtained  by  the  angular  momentum 
matrix  expansion.  The  use  of  the  Wang  representation  of  the  [3  x  (2J  +1)] 
dimensioned  matrices  for  the  Y^  symmetry  (v^  or  v^)  levels  into  four  blocks. 

In  this  fractorization,  three  of  the  four  blocks  have  identical  eigenvalue 
and  are  associated  with  and  F2  symmetry  rotational  levels  whereas  the 
remaining  block  contains  states  of  A^ ,  A^  and  E  symmetry. 


V 
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The  block  form  of  the  rotational  Hamiltonian  coupling  and 
vibrational  levels  is  given  by: 
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The  numerical  results  for  Che  ground  state  of  the  methane  molecule 
are  given  In  Appendix  A.  It  was  found  that  angular  momentum  values  up  to 
J  %  40  could  be  solved  with  modest  amounts  of  computer  time. 

The  Jellium  model  for  high  symmetry  molecules  was  developed  in  order 
to  simplify  the  calculation  of  the  coupled  yibrations-rotation  Hamiltonian. 
This  simplification  is  both  in  terms  of  limiting  the  number  of  parameters 
for  the  Hamiltonian  and  in  the  diagonalization  of  the  resulting  matrix. 

These  questions  together  with  a  complete  exposition  of  the  model  are  the 
subject  of  a  paper  which  is  included  as  Appendix  B.  This  paper  was  submitted 
to  the  Journal  of  Molecular  Spectroscopy. 

Finally,  during  the  final  few  months  of  this  contract,  consideration 
was  given  to  a  calculation  of  the  vibrational  frequencies  of  the  methane 
molecule.  A  force  constant  model  for  methane  was  developed  and  the  evaluation 
of  the  force  constants  based  on  the  electronic  structure  was  considered. 

Our  electronic  state  calculations  are  carried  out  using  scattered  wave  X 

a 

programs  developed  at  MIT.  This  work  was  carried  out  in  collaboration  with 
Miss  Iris  Howard  who  is  the  holder  of  an  IBM  Graduate  Fellowship.  Work  is 
continuing  on  this  project  under  an  AIUSR  grant. 
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Appendix  A.  Ground  State  Rotational  Levels  for  Methane 

The  rotational  Hamiltonian  for  the  methane  molecule  with  symmetry  is 
diagonalized  in  a  fortran  program  which  is  available  on  request.  The  program 
uses  as  inputs  the  coefficients  of  the  terms  in  the  Hamiltonian  and  computes 
the  eigenvalues  for  various  values  of  angular  momentum,  J.  The  symmetry  of 
the  eigenvector  is  also  determined. 

The  program  makes  use  of  the  Wang  transformation  on  the  angular  momentum 
matrices  of  the  rotational  Hamiltonian  into  four  blocks.  Because  of  the 
symmetry  imposed  degeneracies,  only  two  of  these  blocks  need  to  be  diagonalized. 

The  output  from  this  computer  program  for  the  ground  state  of  the  methan. 


molecule  is  given  below. 
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Appendix  B.  Jellium  Model 

The  jelliura  model  is  developed  theoretically  in  the  preprint  which  is 
included  below.  A  computer  program  based  on  this  model  was  also  written.  Hov 
ever,  the  program  used  to  test  the  model  for  the  coupled  v^»  bands  in 
methane,  as  well  as  the  experimental  data  set  used  to  determine  the  goodness 
of  fit  was  that  of  Professor  Alan  Robiettle  of  the  Department  of  Chemistry, 
Reading  University,  Reading,  England. 
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Abstract 


Group  theoretical  techniques  are  employed  to  derive  the  form 
of  the  most  general  phenomenological  rotational  Hamiltonian  associated 
with  a  given  molecular  vibrational  state  or  set  of  states.  The 
expansion  is  given  in  terms  of  the  symmetrized  products  of  angular 
momentum  matrices.  A  perturbation  expansion  for  the  rotational 
states  is  developed  based  on  tetrahedral  distortions  of  a  spherically 
symmetric  distribution  of  mass.  This  jellium  approximation 
significantly  reduces  the  number  of  parameters  in  the  rotational 
Hamiltonian.  The  specific  application  to  methane  is  indicated. 
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I.  Introduction 

Phenomenological  Hamiltonians  are  widely  used  in  physics  to  model 

physical  systems.  Their  introduction  is  motivated  by  either  (1)  an  attempt 

to  model  an  extremely  complex  many-particle  system  by  a  simple  one-particle 

mathematical  formulation  or  (2)  the  necessity  of  refining  theoretical  models 

to  fit  a  multitude  of  precise  experimental  data.  An  example  of  the  first 

point  is  the  introduction  of  pseudopotentials  in  the  theory  of  the  electronic 

structure  of  solids ,  ^  whereas  an  example  of  the  second  is  the  use  of  phenomeno- 

2 

logical  rotational  Hamiltonians  to  fit  the  IR  and  Raman  spectra  of  molecules. 

The  development  of  these  Hamiltonians  has  traditionally  made  extensive 

use  of  symmetry  in  determining  the  functional  form.  The  high  symmetry  molecul 

require  extensive  use  of  group  theory  for  the  primary  reason  that  symmetry 

imposes  mode  degeneracies  which  must  be  explicitly  included  in  order  to 

successfully  model  the  vibrational  and  rotational  states.  The  extensively 

3 

studied  molecule  with  T^  symmetry,  CH^,  will  be  explicitly  considered  in 

this  paper.  However  the  higher  symmetry  molecules  with  0^  symmetry,  e.g. 

SF^  or  UF^  can  be  treated  by  exactly  the  same  formalism  as  that  developed  here 

Tre  methane  molecule  has  4  distinct  vibrational  frequencies  labelled  by 

2 

Vl’  V2’  V3  8nC*  V4’  corresPon<^in8  t0  the  9  degrees  of  vibrational  motion, 
which  consist  of  one  non-degenerate  mode  with  A^  symmetry,  two  degenerate 
modes  with  E  symmetry  and  two  three-fold  and  modes  with  F2  symmetry. 
In  ad'ltion  to  these  symmetry- imposed  degeneracies  (which  require  special 
consideration  in  the  derivation  of  the  phenomenological  Hamiltonian),  there 
are  near-degeneracies  that  are  accidental.  Two  examples  of  such  near¬ 
degeneracies  in  methane  are:  (1)  V2  and  v^,  and  (2)  v^,  the  overtones 
(2V2>  and  (2v^)  and  the  combination  modes  (V2  +  v^).  The  fitting  of  the 
vibronic  modes  in  CH^  is  greatly  perturbed  by  the  interactions  between  the 
levels  of  these  degenerate  and  nearly  degenerate  states. 


3 


This  paper  addresses  the  formulation  of  the  phenomenological  Hamiltonian 
for  the  rotational-vibrational  states  in  a  high  symmetry  molecule  and  develops 
a  perturbation  expansion  based  on  deviations  from  spherical  symmetry.  This 
formulation  makes  it  possible  to  treat  a  collection  of  degenerate  and  nearly 
degenerate  levels  using  the  minimal  number  of  terms  required  by  symmetry 
considerations  in  the  phenomenological  Hamiltonian.  Using  the  functional  form 
based  on  group  theoretical  considerations,  the  coefficients  of  the  model 
Hamiltonians  can  be  evaluated  by  comparison  with  experiment.  Though  application 
is  made  to  CH^  having  modes  with  A^,  E  and  2?^  symmetries,  the  same  formulation 
is  easily  extended  to  other  molecules  with  T^  symmetry,  such  as  Ni(CO)^  having 
modes  with  2A^,  2E,  and  4F£  symmetries.  The  extension  to  high  symmetry 
molecules  with  different  cubic  symmetries  can  also  be  carried  out  based  on  the 
functional  form  presented  in  this  paper. 

The  general  symmetry  considerations  which  apply  to  all  phenomenological 
Hamiltonians  are  discussed  in  Section  II  and  these  principles  are  applied  to 
vibrational- rotational  Hamiltonians  for  Td  symmetry  molecules.  Section  III 
develops  a  perturbation  expansion  based  on  the  spherically  symmetric 
phenomenological  Hamiltonian  for  "jellium"  which  allows  one  to  greatly  reduce 
the  number  of  expansion  coefficients.  Section  IV  makes  the  application  of  the 
perturbed  jellium  model  to  the  methane  molecule  and  finally  Section  V  gives 
a  brief  summary  of  the  current  paper. 
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II .  General  Theoretical  Development 
A  phenomenological  Hamiltonian 


ft 


ri 


(0  ,  0  ,  ...  0  )  in  general  depends 


on  the  operators  0  and  is  a  scalar  with  respect  to  the  symmetry  operations 
of  the  symmetry  group  of  the  system.  Thus  if  R  is  any  symmetry  operation 
(including  time  inversion),  that  leaves  the  Hamlitonian  invariant;  then 

"Kr/V  e2-  -V  ‘  tfr/V  VV-  (1) 

The  invariance  of^  under  the  symmetry  operations  R  does  not  otherwise 
restrict  the  form  of  $rv  so  that^f  can  be  a  matrix,  which  is  the  case 


if 


ft. 


represents  the  rotations  of  a  vibrating  molecule.  The  operators  0. 


are  also  not  necessarily  scalars  under  the  symmetry  operations  R.  These 

apparent  complexities  can  be  systematically  handled  by  elementary  group 

theoretical  considerations,  and  this  is  the  focus  of  the  present  paper. 

The  symmetrized  form  of  Arp,  (0,,  0„,  ...6  )  can  be  based  on  the 

v  1  l  l  n 

4 

assumption  that  its  functional  form  is  a  power  series  in  one  or  more  operators. 
This  assumption  is  not  necessary,  but  is  convenient  because  a  function  of  an 
operator  f(0)'iJ  is  usually  defined  in  terms  of  a  power  series  expansion  in  the 
operator,  i.e.  (f(0)  +  f'(O)0  From  this  point  of  view,  the  Hamiltonian 

is  written  in  terms  of  a  complete  set  of  basis  matrices  B(I\)  which  transform 
as  the  irreducible  representation  T, 


r  (0  ,  0  , ...6  )  =  I{B(r),  Ar  (6.  .... ,0  )}r 
‘1  i  l  n  i  l^l  n*l 


(2) 


where  the  Ap^(0^...0n)  are  symmetrized  polynomial  functions  of  the  operators 
0^,  02»...0n  that  transform  according  to  the  irreducible  representation 
The  basis  matrices  B(T^)  can  be  classified  according  to  the  irreducible 
representation  of  the  symmetry  group  of  the  molecule  and  the  scalar  combination 
of  B  and  A  is  taken  in  Eq.  (2).  Time  reversal  symmetry  must  also  be  explicitly 
considered,  so  that  the  basis  matrices  are  either  even  or  odd  under  time  reversal. 
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If  a  given  irreducible  representation  ia  n-dimensional,  then  there  are  n2 
basis  matrices  of  rank  n,  of  which  n(n-l)/2  are  pdd  under  time  reversal  and 
n(n+l)/2  are  even.  The  formulation  given  here  is  quite  general.  When  applied 
to  the  development  of  phenomenological  Hamiltonians  for  solids,  space  groups 
must  be  used,  whereas  for  molecules  point  groups  are  used. 

A  general  expansion  for  degenerate  or  nearly-degenerate  states  (whether 
vibrational  or  electronic)  is  obtained  by  use  of  angular  momentum  matrices?  L. 
Use  of  these  basis  matrices  is  common  for  electronic5 E  vs  k  or  phonon6w  vs  q 
in  crystals  and  was  first  introduced  for  molecular  rotational  Hamiltonians  by 
Landau  and  Lifshitz.7  The  rotational  Hamiltonian  for  the  vibrational  state 
characterized  by  the  angular  momentum  quantum  number  l  taken  to  order  (2n)  in 
the  angular  momentum  matrices  J  is  designated  by  K  <2n)  and  is  given  by 


2n  U  ,  } 


^(2n)(j)  -  ll  l  l  a£U)  (u,X,v)[j(j+l)]V{jfX)(L(ll)),  j£W"2v)'^ 

jj-o  v»o  * 


(A£) 

The  parameters  of  the  expansion  ap^  (y,X,v)  are  indexed  according  to  symmetry 

type  T±  and  to  the  three  indices  y,X,  and  V,  where  y  denotes  the  exponent  of 

angular  momentum  J,X  the  exponent  for  angular  momentum  \  and  v  the 

exponent  for  the  product  (J  •  J)»  J (j+1)  where  j(j+l)  is  the  eigenvalue  of  J  . 

(X')  -t 

The  symmetrized  combinations  of  angular  momentum  matrices  Jp^  (J)  are 
labelled  according  to  symmetry  type  and  to  the  order  of  the  basis  function  X'. 
Explicit  symmetrized  linear  combinations  corresponding  to  T^  symmetry  are  given  by 
Altmann  and  Cracknell®  for  0$X'<12.  The  curly  bracket  (  ,  indicates  that 
symmetrized  combinations  of  Jp  ^(L^)  and  jp^^”  2V\5)  are  taken  so  that 
a  scalar  Hamiltonian  results  and  the  non-commuting  quantum  mechanical  operators 
are  symmetrized.  The  prime  on  the  summation  over  X  denotes  the  constraint 
that  X  +  y  •  even  and  the  sum  over  v  is  cut  off  at  the  largest  integer  in  (y/2). 


Further  constraints  on  the  summations  in  Eq.  (3)  are  presented  below.  These 

.  (.A 


constraints  depend  on  the  properties  of  the  symmetrized  basis  functions  jp 
which  are  now  discussed. 


The  symmetrized  products  of  angular  momentum  matrices  ^ )  span 


the  vector  space  appropriate  to  angular  momentum  £  and  are  determined  by  the 
dimension  of  the  set  of  levels.  For  example,  for  a  non-degenerate  level,  only 
one  matrix  is  needed,  the  unit  matrix  (1),  so  that  =  (1) 

for  X  =  0.  For  a  3-fold  degenerate  level,  l  =  1  in  Eq.  (3)  and  9 
basis  functions  are  required  for  span  the  vector  space.  These  basis  functions 


are  found  from  the  direct  product  of  the  symmetry  type  F^  for  l-l,  namely 


x  F^  =  +  F^  +  E  +  F2*  These  9  basis  functions  for  T^  symmetry 

are  given  by  where  SL  «=  1  and  X  =  0,1,2.  The  matrices 

given  by^* 


(1) 


(1) 


(1) 


respectively,  generate  the  9  basis  functions  for  a  3-fold  degenerate  level 
with  F,  symmetry  in  the  T^  group.  The  normalization  is  given  by 


l  <j‘0<J>,  jfX)  (J))  .  =  <5  •  J)X  (5) 

r  i  i  1 

i 

where  the  sum  is  over  ail  the  irreducible  representations  that  occur  for 

2  _ 

a  given  a,  J  #J  =  J  and  the  curly  bracket  denotes  a  scalar  product  of  the 

symmetrized  basis  functions.  For  X  £  5,  multiple  occurrences  of  a  given 
irreducible  representation  are  found,  as  for  example  F^  and  F^,  for  X  *  5, 
and  F^  and  F2,  for  X  «  6.  Thus,  if  X  or  (y-  2v)  in  Eq.  (3)  exceed  A,  then 
the  sum  over  representations  must  include  multiple  occurrences  such  as  F^ 
and  Fj,  for  X  ■  5,  and  F^  and  F^,  for  X  -  6.  Similarly  if  l  >  2  in  Eq.  (3), 
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then  X  will  exceed  4  and  multiple  occurrences  of  representations  must  again 
by  explicitly  included.  Also  note  that  the  functions  of  (l/^)  for  even 

X  are  even  under  time  reversal,  while  the  functions  with  odd  X  are  odd  under 
time  inversion. 

Application  of  the  general  formulation  is  made  to  the  various  common 
syranetry  types  for  vibrational-rotational  levels  occurring  for  molecules  with 
symmetry.  In  particular  for  1=0  (representing  the  ground  state  or  a  non¬ 
degenerate  vibrational  state  with  either  A^  or  symmetry  (i,e,  the  v^ 
vibrational  state  of  methane  which  has  A^  symmetry),  then  the  only  basis 
matrix  which  enters  Eq.  (3)  is  =  1.  Thus  the  sum  over  symmetries 

is  restricted  to  those  combinations  of  (J)  having  A^  symmetry.  Time 

inversion  invariance  further  restricts  the  sum  to  p  =  even  integer  which  is 

written  as  2m.  Thus  the  most  general  rotation  Hamiltonian  for  a  non-degenerate 
9,10 


level  is 


K 


(2n)  __  n  m 

<J>  "  l  l 


a<°0)  (2m,0,v)  [j(j+l)]V  jA  (2(m"V))(J). 
m=0  V=0  A1  1 


(6) 


Table  1  gives  the  number  of  aA^  coefficients  which  occur  in  Eq.  (3) 

for  1  values  in  the  range  0  £  1  £  3,  0  £  V  £  10  and  all  allowed  values  of  X,v. 

In  particular,  for  1=0,  then  X  is  restricted  to  X=0,  and  V  is  restricted  to 

0  £  V  i  m.  Table  1  then  gives  one  non-vanishing  term  for  m=0  (namely  with 

coefficient  aA^  (0,0,0),  and  one  non-vanishing  term  for  m=l  (namely  with 

coefficient  a  A^°^ (2 *0, 1) ) ,  since  the  symmetry  type  A^  does  not  occur  for  X=2. 

Table  1  shows  that  for  m =2  there  are  two  non-vanishing  terms  in  Eq.  (6),  and 

these  correspond  to  coefficients  a^®’0^ (4,0,0)  and  a  A^0,0^ (4,0,2) .  The 

rotational  Hamiltonian  "5^  00°^ (^>  obtained  from  E9*  (6)  has  precisely  the 

9  10 

same  form  as  the  Hamiltonian  used  previously  by  other  workers,  ’ 


The  next  important  application  of  Eq.  (3)  to  the  vibrational-rotational 

spectrum  for  methane  is  to  the  and  levels  which  have  F2  symmetry^  and 

are  3-fold  degenerate  ($,=1).  The  Hamiltonian  for  £=1  can  be  written  for 

11 

levels  with  either  or  F2  symmetry  as: 

0  V 11  m*=0  V=0  A1  1  1  '  1 

n-l  tn  (2 (m-v)+l) 

+  11  a  (11)(2nH-l,l,V)[j(j+l)]  {j*1'  (L(1)),  j  (3)  } 

m=0  V=0  F1  F1  F1  A1 

n  m_1  nn  „  mn  _*/ (2 (m-v) ) 

+  l  l  a(F1)(2m,2,v)[j(j+l)]V  (j^2)(L(1)),  3  (3)  } 

m=l  V=0  E  EE  Aj 


l  l  a  i.11)  <2m,2,v)[j(j+l)]V  { j <2) (L(1) ) ,  (3)  ^ )  .  (7) 


m=l  v=0  2 


in  which  the  sum  on  X  in  Eq.  (3)  has  been  written  explicitly.  The  first  sum, 

corresponding  to  A=0, is  identical  to  the  Hamiltonian  in  Eq.  (6),  except  that 

JA]^(l^)  is  a  (3x3)  unit  matrix,  while  Ja^(L^)  is  the  scalar  1.  The 

second  term  in  Eq.  (7)  is  for  X-l,  and  time  reversal  symmetry  restricts  the 

sum  to  odd  powers  of  J  as  indicated  since  only  T^  =  F^  symmetry 

occurs  for  jp^^(L^).  To  obtain  a  scalar  product  with  A^  symmetry,  the 

basis  function  j  2V^(J)  must  have  T  =  F  .  The  number  of  coefficients 

occurring  for  each  order  m  is  found  in  Table  1  for  £=1  and  symmetry. 

Similary  for  '’2,  two  symmetry  types  occur,  namely  E  and  F^  and  time  reversal 

symmetry  requires  even  powers  of  J  in  the  last  two  sums  in  Eq.  (7).  The 

number  of  expansion  coefficients  is  again  found  from  Table  1  using  the 

entries  listed  under  £*1,  and  E  and  F2  symmetries,  respectively.  The 

11  12 

Hamiltonian  in  Eq.  (7)  is  identical  to  that  used  by  Robiette  et  al.  ’ 
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To  treat  the  \>2  level  in  methane,  it  is  necessary  to  write  Eq.  (3)  for 
which  is  appropriate  for  treating  the  coupled  E  and  F2  vibrational  modes.  An 
isolated  E  mode  can  be  treated  by  the  introduction  of  the  4  basis  matrices  jp^^(d  ) 
having  ExE  *  +  A^  +  E  symmetry  and  given  by 


.(0),.  v  ,1  0  ,(1),.  .  .1  0. 

\  <de}  (0  1}*  jA2  (de)  =  (0  -1>  ’ 


,  (2),.  ,  ,0  1.  .  .(2),.  .  ,0  0. 

^de^  “  ^0  O’  and  ^E*  ^de^  “  ^1  0^ 


(8) 


where  the  8  index  is  used  as  a  convenient  label  to  indicate  time  inversion 

and  dft  refers  to  the  part  of  the  d-f unction  that  has  E  symmetry.  The  resulting 

Hamiltonian  appropriate  to  the  v2  vibrational  state  of  methane  can  then  be 
11 

written  as 


<2n>  -  «  “  (ee-)  v  (2(m-v)) 

“  I  I  a  Ce)  (2m,0,v)  {j  (j-fl)  ]V  {jA(0)(d  ),  j.  (J)  } 


ee 


m=0  v=0  1 


'Ai  '~e' 


n  m 

+  l  l 


(ee)  v  fl)  (2(f-V)+l) 

alee;(2m+l,l,v)[j(j+l)]V{jy)(d  ),  J  (j)  } 

1  v=0  A„  2  2  A1 


m 


+  l  I 

m=l  v=0 


*2 

a(ee) 

E 


V  (2)  (2(m-v) ) 

(2m,2,v) [j (j+1) ]V( jg  ' (dg>  >  JE  tf)  }  A 


(9) 


in  which  the  proper  linear  combination  of  operators  J  has  been 
taken  to  preserve  time  inversion  symmetry.  The  number  of  coefficients 

ap£ee^(p, 8, v)  corresponding  to  each  value  of  V  is  listed  in  Table  2  under 

r  *  E.  r  -  E  and  T  -  A. ,  A„  or  E.  Further  discussion  of  this  table  is 
i  ’  j  k  1  2 

given  below.  The  Hamiltonian  given  by  Eq.  (9)  has  also  been  used  previously 
to  treat  the  v^  level  of  methane. 

The  present  work  makes  it  possible  to  simultaneously  derive  the  coupling 
between  degenerate  and  nearly  degenerate  vibrational  levels  to  various  orders 
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in  J.  The  general  formulation  for  this  coupling  is  given  below  in  the  context 
of  the  coupling  between  an  E  and  an  F£  (or  F^)  vibrational  level.  This 
particular  coupling  block  is  important  for  the  methane  molecule  and  is  present 
in  the  Hamiltonian  ^^2°^  ^  for  contains  the  block 


given  by  Eq.  (9).  In  Section  in  is  presented  a  perturbation  approach  that 
reduces  the  number  of  expansion  coefficients  that  must  be  used  to  yield  energy 
eigenvalues  for  a  given  order  in  the  expansion. 

The  full  expansion  for  1=2  in  symmetry  reduces  to  coupled  E  and  F^ 

(or  E  and  F^)  vibrational  states.  The  Hamiltonian  for  £=2  is  then  written  as: 


11 


# 


tee"'*  I 


(2n) 

22 


(J)  = 


files'  #«*> 


(J) 


(10) 


where  the  diagonal  blocks  are  given  by  Eq.  (9)  for /J^ge  and  Eq .  (9)  for 


V. 


The  diagonal  blocks  yield  energies  for  uncoupled  v.  and  v, 

2  “4 


levels  for  methane.  Since  these  levels  are  nearly  degenerate,  they  are 
strongly  coupled,  and  this  coupling  is  treated  by  the  off-diagonal  blocks 


« 


^n^(J).  The  symmetry  and  functional  form  for  these  coupling  blocks  are 


found  by  expanding  the  general  Hamiltonian  in  Eq.  (3)  and  then  identifying 

the  resulting  off-diagonal  blocks  with^^n^  (J)  .  In  this  expansion  the 
-*•(2) 

basis  matrices  L  are  5x5  matrices  appropriate  to  l  =  2.  The  number  of 
expansion  coefficients  for  each  value  of  y  (i.e.  Jw)  is  found  in  Table  2 
under  ■  E,  T  «  F2.  In  this  table  blocks  with  the  same  symmetry  are  coupled 

O  2 

foru-  0  (i.e.  J  )  and  these  terms  are  called  Fermi  coupling  terms.  the  coupling 


of  two  different  blocks  by  terms  linear  in  J  are  called  Coriolis  coupling 
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terms. ^  Terms  in  the  coupling  block  (J)  corresponding  to  y  5  2  have  no 

special  designations,  as  do  the  terms  corresponding  to  y  =  0,1.  From  Table  2 
it  is  seen  that  Eq.  (10)  has  no  Fermi  coupling  terms  and  1  Coriolis  coupling 
term.  Table  2  also  gives  the  number  of  terms  which  occur  in  the  coupling 
blocks  for  each  y  in  the  range  2  s'  y  £  10.  Also  included  in  Table  2  is  the 
number  of  expansion  coefficients  associated  with  general  coupling  blocks  of 
T±  and  symmetries;  such  coupling  blocks  occur  in  Hamiltonians  Kr  J) 
for  £  >  2. 

The  rotational  Hamiltonians  and  (3)  presented 

in  this  section  have  all  been  used  to  fit  the  rotational  levels  for  the 

vibration  states  in  the  methane  molecule.  The  total  number  of  coefficients 

which  is  required  for  the  2nth  order  phenomenological  Hamiltonian  is  given 

in  the  first  3  lines  of  Table  3.  Thus  we  see  that  for  the  uncoupled  v^,  V£, 

Vj,  and  levels  the  8th  order  Hamiltonians  have  11,26,48  and  48  expansion 

coefficients  respectively.  The  applications  of  phenomenological  Hamiltonians 

that  have  previously  been  made  to  methane  have  not  been  able  to  fit  the 

experimentally  measured  levels  over  the  full  range  of  j  quantum  numbers 

12 

that  are  experimentally  accessible.  Various  schemes  have  been  proposed  to 
improve  the  agreement  between  the  experimental  data  and  the  phenomenological 
Hamiltonians,  and  these  schemes  have  generally  involved  inclusion  of  coupling 
terms  between  nearly  degenerate  vibrational  states.  ^he  number  of  coefficients 
listed  for  the  coupled  vibrational  levels  on  the  last  3  lines  of  Table  3 
assumes  that  the  coupling  terms  are  included  to  the  same  order  in  J  as  in  the 
diagonal  blocks.  In  application  of  these  equations  to  a  physical  problem, 
the  magnitudes  of  the  coupling  parameters  are  such  that  they  often  can  be 
treated  in  perturbation  theory  following  the  formulation  given  here.  In 
addition,  a  somewhat  different  type  of  perturbation  theory,  based  on 


approximate  spherical  symmetry  for  the  molecule  (a  jellium  approximation  for 
the  unperturbed  molecule),  can  be  useful  in  certain  physical  situations.  This 
type  of  perturbation  approach  is  presented  in  the  next  section. 

III.  Jellium  Perturbation  Model 

The  formulation  of  the  phenomenological  Hamiltonian  in  Section  II, 

though  mathematically  correct,  is  often  inconvenient  to  apply  because  of 

CZl) 

the  large  number  of  coefficients  ap^  (y,\,v)  that  must  be  determined  from 
experiment.  In  this  section,  a  perturbation  theory  approach  is  presented 
which  yields  a  more  convenient  form  of  the  phenomenological  Hamiltonian 
specific  to  a  given  molecule.  In  this  formulation,  a  jellium  model  with 
spherical  symmetry  is  assumed  for  the  zeroth  order  Hamiltonian,  and  the 
specific  symmetry  is  introduced  by  a  perturbation  term  with  the  proper 
symmetry.  For  example,  the  perturbation  theory  for  a  molecule  with  Td 
symmetry  would  be  written  as 


(11) 


in  which  the  unperturbed  zeroth  order  Hamiltonian  has  full  rotational 

symmetry  and  is  written  as  and  the  perturbation  Hamiltonian  is 

written  as  A  discussion  of  each  term  in  Eq.  (11)  is  now  presented. 

The  assumption  of  spherical  symmetry  for  imPlies  that  a11 

SLi. 

expansion  parameters  (y,A,v)  in  Eq.  (3)  for  a  given  order  y  are  equal,  independ 
of  symmetry  type  so  that  the  expansion  parameters  for  the  unperturbed 
term  are  written  as 


lar**£)  (M.*,\>)1  =  6 


A,y-2v 


0  <W)(V) 


This  assumption  and  Eq.  (3)  enables  one  to  write  the  unperturbed  Hamiltonian  as 
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[^un)(3)]0  “  r  {yTlo>(U)(v)  lj(j+l)]-v  (px<L(t> 


J) 


(13) 


in  which 


X=0  v=0 

J)  is  defined  by 


(L 


sv 


J)  - 1 

r. 


<4A)  tw),  4X)  (d))A 

ri  ri  ai 


(14) 


Using  the  properties  of  spherical  harmonics,  the  sum  in  Eq.  (14)  can  be 
carried  out  to  give 


[*a+i)j(j+i)]A/2 


( 


KA(£+l)j(j+l)] 


1/2 


j  (15) 


where  P^(z)  is  the  Legendre  polynomial  of  order  X.  We  further  note  that 

when  =  J  in  Eq.  (15),  the  function •  J)  reduces  to  (J  •  J)\ 

consistent  with  Eq.  (5).  The  jellium  approximation  thus  leads  to  level 

orderings  and  degeneracies  appropriate  to  vector  sums  of  the  angular  momenta 

Lv  +  J  with  quantum  numbers  j+&,  j  +S.-1, . .  .  j-&.  The  (2£+l)  x  (2)1+1) 

Hamiltonian  defined  by  Eq.  (13)  has  precisely  the  same  block  form  as  the 

Hamiltonian  of  Eq.  (3)  with  T^  symmetry,  except  that  the  expansion  coefficients 

zt 

ap^(p,X,v)  are  constrained  by  Eq.  (12).  In  this  way,  a  major  reduction  in 

the  number  of  expansion  coefficients  is  achieved.  The  number  of  independent 

(.It) 

coefficients  (v)  in  Eq.  (13)  is  listed  in  Table  4  according  to  index  i 

and  order  of  the  Hamiltonian  2n,  where  X  $  2t  and  v  is  the  largest  integer 

less  than  or  equal  to  (n-X/2).  As  an  example,  for  t=0,  truncation  of  the 

Hamiltonian  at  results  in  four  expansion  parameters,  written  as  o  ^^^(v); 

o 

v**0,l,2,3.  To  obtain  off-diagonal  blocks  in  the  Hamiltonian  with  all  the 
distinct  symmetry  types  that  occur  in  T^  symmetry  when  the  perturbation 
term  in  Eq,  (11)  is  introduced,  the  expansion  of  the  jellium  Hamiltonian  of 
Eq,  (13)  must  be  carried  to  terms  n  s  *. 
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Referring  again  to  Eq.  (11),  the  perturbation  Hamiltonian  ( 

can  be  written  explicitly  using  Eq.  (3)  for  0  <  p  <  n.  To  yield  contribution: 

to  the  energy  of  order  2n,  the  perturbation  Hamiltonian  is  truncated  at  order  n 

Thus,  when  the  perturbation  energy  is  calculated  for  levels  that  are  split  by 

an  interactions  with  symmetry,  the  coupling  blocks,  taken  to  order  Jn  on 

the  off-diagonal  positions  of  the  Hamiltonian,  yield  contributions  to  the 

2n 

block-diagonal  terms  of  order  J  ,  which  is  of  the  same  order  as  the  jellium 

Hamiltonian.  The  coefficients  in  the  perturbation  Hamiltonian  are  written  as 

( n) 

aFl  (u;X,v),  and  are  cataloged  by  symmetry  type  r  ,  in  the  same  way  as  the 

(iZ)  (ZZ) 

ap^  (u,X,v)  coefficients  of  Eq.  (3).  The  ap^  coefficient  however  dffers 

(£.£.)  (ZZ) 

from  the  ap^  coefficients  insofar  as  the  ap^  are  constrained  to  yield 

(i  Z) 

the  spherically  symmetric  a '  (v)  when  summed  over  all  symmetry  types 

l  ar(U)(A  +  2v,  X,  v)  «=  n  ax<M)(v)  (16) 

r  Ai  r  A 

‘i 

where  is  the  number  of  distinct  representations. 

The  large  reduction  in  the  number  of  expansion  coefficients  for  the 

coupled  vibrations  represented  by  Eq.  (11)  relative  to  the  general  form  of 

Eq.  (3)  is  shown  explicitly  in  Table  5.  In  particular  the  Jl=2  Hamiltonian 

appropriate  to  the  coupled  Vj,  vibrational  states  in  methane  can  be 

modelled  using  Eq.  (11)  with  39  coefficients  to  8th  order  which  is  to  be 

contrasted  to  114  coefficients  for  the  full  T,  symmetry  expansion  to  8th 

a 

The  jellium  perturbation  theory  developed  in  this  section  gives  a 
major  reduction  of  expansion  coefficients,  especially  for  Z  >  2  as  shown 
in  Table  5,  These  higher  angular  momentum  states  are  directly  applicable 
to  combination  and  overtone  vibrational  states  in  methane  as  is  detailed 
In  the  next  section. 
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IV.  Proposed  Application  of  Jelllum  Model  to  Methane 

Based  on  the  perturbation  theory  developed  In  the  last  section,  we 
now  proceed  to  outline  an  approach  to  the  Hamiltonian  for  the  methane  levels 
which  focusses  on  those  groups  of  states  that  lie  in  approximately  the  same 
frequency  range.  An  application  of  the  jellium  perturbation  expansion  to  the 
coupled  modes  in  methane  was  carried  out  by  Prof.  A.G.  Robiette  using 

his  collected  experimental  data  set  for  these  levels.  The  results  are 
quoted  in  the  footnote  to  the  table  5,  Prof.  Robiette  achieves  somewhat 
better  fits  to  the  experimental  data  by  judicious  selection  of  certain  terms 
in  the  exact  Hamiltonian  of  Eq.  (10). 

All  the  vibrations  of  the  jellium  molecule  can  be  approximately 
specified  in  terms  of  the  three  fundamental  vibrations  which  are  specified 
by  l  =  0.1  and  2.  The  state  of  the  jellium  molecule  is  then  denoted  by  a  set 
of  3  integers  n^  =  (n^,  n^,  n^)  which  specify  how  many  quanta  n^  are  present 
in  each  of  the  three  l  states.  When  either  harmonics  and/or  combination 
modes  are  present,  it  is  necessary  to  use  a  vector  model  to  obtain  the  total 
vibrational  angular  momentum  quantum  number  .  The  jellium  model  Hamiltonian 
is  written  in  terms  of  ,  as  are  the  corresponding  energy  levels  given  on 
the  left  side  of  Table  6  and  denoted  by  n^  and  .  We  note  from  Table  6 
that  in  the  methane  molecule  there  are  a  number  of  accidental  near-degeneracies; 
more  specifically  many  of  the  vibrations  are  such  that  the  methane  levels 
occur  in  approximate  multiples  of  1450  cm  \  as  emphasized  by  the  grouping  of 
levels  given  in  Table  6.  Thus  the  approximate  position  of  the  jellium 
vibrations  in  methane  follow  the  empirical  relations  given  by 


u^Cn^)  :  1450  (n(nQ 


nl> 


n2] 


cm 


-i 


(17) 
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Thus  the  ground  state  (000)  forj^  “  0  is  isolated,  doesn't  interact  strongly 
with  any  other  state  and  remains  non-degenerate  when  the  perturbation  with 
Tj  symmetry  is  introduced.  The  (001)  jellium  state  for ^  =  2  has  5 
degenerate  vibrational  levels  at  1450  cm  ^  and  gives  rise  to  a  coupled 
Hamiltonian  having  the  form  of  Eq.  (13)  and  this  level  splits  under  the 
perturbation  of  the  form  of  Eq.  (1)  into  a  3-fold  level  with 
symmetry  and  a  higher  lying  2- fold 

level  with  E  symmetry.  The  next  set  of  states  at  about  2900  cm  *  consists 
of  the  (100)  level  with<^  =  0,  the  (010)  3-fold  degenerate  level  with^=  1 
and  the  (002)  level  with^£  =  4,2,  and  0.  The  (002)  overtone  for  £  =  2 
would  according  to  the  vector  model  correspond  to  =  4, 3, 2, 1,0  but  the 
odd  values  of  do  not  occur  because  these  states  are  odd  under  the 


permutation  of  identical  particles,  and  phonons  are  Bose  particles.  The 
Hamiltonian  for  the  coupled  jellium  states  with  energies  near  2900  cm  ^ 


is  written  as 


n£,  n£' 


1  O' 2900^0  jHi  X'  |  8) 

13  -*• 

containing  5  x  5  =  25  blocks  which  are  labelled  by  the  n^, ^  indices 
listed  in  Table  5  and  enumerated  explicitly  above.  The  most  important 
coupling  terms  are  those  for  (which  contains  Fermi  coupling  as  the 

leading  term  in  the  expansion)  andj£=  ±  1  (where  the  leading  term  in 
the  expansion  represents  Coriolis  coupling).  The  introduction  of  the 
perturbation  with  T^  symmetry  follows  the  perturbation  theory  of  Eq.  (11). 

The  resulting  splitting  of  the  jellium  levels  on  the  left  hand  columns  of 

Table  b  are  indicated  in  the  right  hand  columns  of  Table  6,  in  terms  of 

the  appropriate  irreducible  representations  of  T^,  the  number  of  phonons  n  ^ 

with  mode  frequencies  vj^*V2,V3  an<*  V4  ®nd  the  experimental  values  of  the  mode 

frequencies^  • 
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The  groupings  of  levels  near  4350  cm  ^  are  also  indicated  in  Table  6, 
together  with  the  jellium  model  quantum  numbers  'and  mode  frequencies 
«  (n^) <  and  also  the  corresponding  level  splittings  which  occur  in 

symmetry. 

V.  Summary 

High  symmetry  molecules  represent  a  theoretical  challenge  in  molecular 
physics  for  two  essential  reasons:  (1)  Molecules  with  cubic  symmetry  must 
have  at  least  4  identical  atoms  (perhaps  more)  and  hence  6  or  more  internal 
vibrational  degrees  of  freedom  (i.e.  a  large  number  of  modes)  and  (2)  symmetry 
forces  some  of  these  nodes  to  be  degenerate.  The  high  symmetry  furthermore 
increases  the  probability  for  near-degeneracies  which  further  complicate  the 
calculation  of  the  vibrational  energies.  The  models  for  these  molecules 
must  therefore  make  extensive  use  of  group  theoretical  techniques. 

This  paper  gives  the  most  general  group  theoretical  form  for  the 
rotational  Hamiltonian  associated  with  the  vibrational  states  of  these  high 
symmetry  molecules  and  the  use  of  this  Hamiltonian  is  illustrated  for  the 
methane  model.  In  addition  a  perturbation  theory  is  developed  based  on  a 
jellium  model  for  themolecule  which  enables  one  to  significantly  reduce  the 
number  of  expansion  coefficients.  The  overtone  and  combination  modes  in 
the  high  symmetry  molecules  occur  in  the  same  frequency  region  as  the 
fundamental  vibrations  and  perturbation  theory  developed  in  this  paper  is 
formulated  to  treat  to  the  coupling  between  this  extended  set  of  nearly 
degenerate  modes  in  the  methane  molecule. 
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Table  2  Number  of  expansion  coefficients  associated  with  the  coupling  terms 
in  the  phenomenological  Hamiltonian  (a,b) 


Number  of  Coefficients 


ri 

Symmetry  y  -  0  1  2 

3 

4 

5 

6 

7 

8 

9 

10 

A1 

A2 

Ai7  Ai 

'  2 

3 

4 

1° 

5 

A1 

A2  A2 

1 

1 

1 

2 

1 

3 

2 

A1 

*2 

El  E 

A 

2 

1 

3 

2 

5 

3 

7 

A1 

A2 

:;]  ’■ 

1 

2 

2 

4 

4 

6 

6 

9 

A1 

A2 

a Fi 

2 

1 

4 

2 

6 

4 

9 

6 

E 

E  Ax  1  1 

2 

3 

4 

l6 

5 

A2 

1 

1 

1° 

2 

1° 

3 

2° 

E  1 

2 

1° 

3 

2° 

5 

3° 

7 

E 

Fi7  Fi  1 

2 

1 

4 

2 

6 

4 

9 

6 

E 

F2i  F2  1 

1 

2 

2 

4 

4 

6 

6 

9 

F1 

Fi?  Ai  1  1 

2 

3 

4 

1° 

5 

F2 

f2J  e  1 

2 

1° 

3 

2° 

5 

3° 

7 

F1 

2 

1° 

4 

2° 

6 

4° 

9 

6° 

F2  1 

1° 

2 

2° 

4 

4° 

4 

6° 

7 

F1 

F2?  A2 

1 

1 

1 

2 

1 

3 

2 

F2 

Fl5  E  1 

2 

1 

3 

2 

5 

3 

7 

F1  1 

2 

1 

4 

2 

6 

4 

9 

6 

F2  1 

1 

2 

2 

4 

4 

6 

6 

9 

a) 

rt  and  denote  the  symmetries  of 

the 

coupled 

terms. 

b) 


The  terms  which  are  odd  under 
zero  superscript  when 

of  the  Hamiltonian,  i.e.  Eqs. 
are  odd  under  time  inversion, 
couples  two  modes  of  the  same 

in  methane)  would  include  the 


the  time  inversion  symmetry  are  labelled  with  a 
These  terms  are  not  allowed  on  the  diagonal  block 

(6),  (7)  and  (9)  do  not  contain  any  terms  that 
However,  the  off-diagonal  interaction  that 
symmetry  (e.g.  the  and  vibrational  levels 

odd  terms  indicated  in  this  table. 


Table  3 


Total  Number  of  Independent  Expansion  Coefficients 

(Hi) 

a_  (p,A,v)  in  the  Phenomenological  T.  Rotational 
1  i 

Hamiltonian  Summed  to  v  *  2n 


Number  of  Independent  Constants 


Vibrational  Symmetry 

Dimension 

2n  =  2 

4 

6 

8 

10 

A1  °r  A2 

lxl 

2 

4 

7 

11 

16 

E 

2x2 

3 

8 

15 

26 

41 

Fx  or  Fj* 

3x3 

5 

13 

27 

48 

78 

(ArE)  or  (A2,E) 

3x3 

6 

15 

29 

51 

81 

(AlfF2)  or  (A2,F1)“? 
(A1,F1)  or  (A2,F2)J 

4x4 

8 

21 

44 

79 

129 

(E,F2)  or  (E,F1) 

5x5 

10 

2S 

62 

114 

189 

a)  The  2nt*1  order  rotational  Hamiltonian  used  by  various  authors  for  the 
Vj  and  states  in  methane  includes  in  addition  the  symmetry 
combinations  with  p  =  2n  +  2.  Hence  their  2, 4, 6, 8  and  10th  order 
Hamiltonians  have  7,16,31,53  and  85  parameters  respectively. 


Table  4 


i  4 

Total 

Number 

of 

Independent  Coefficients  c 

in  the 

Jelllum 

Hamiltonian 

a)  of 

• 

Order 

2n 

Total 

Number 

of 

Coefficients  for 

Order 

2n 

£ 

2n  = 

2 

4 

6 

8 

10 

0 

2 

3 

4 

5 

6 

1 

4 

7 

10 

13 

16 

2 

4 

9 

14 

19 

24 

3 

4 

9 

16 

23 

30 

The  Integers  X  and  V  are  constrained  by  0  6  X  <  2£ 
and  0  <  v  <  {n  in  which  {n  -  X/2)  is  the 

largest  integer  less  than  or  equal  to  n  -  X/2. 


Tab le  5 


Comparison  Between  the  Total  Number  of  Symmetry  Allowed 
Coefficients  in  the  "Full"  Phenomenological  Hamiltonian 
(Eq.  (3))  and  in  the  "Jellium"  Perturbation  Expansion  of 


Eq. 

(11). 

Total 

Number  of 

Coefficients  for 

Order  2n 

J l 

Type 

3n  -  2 

4 

6 

8 

10 

0 

full 

2 

4 

7 

11 

16 

jellium 

2 

3 

4 

6 

7 

1 

full 

5 

13 

27 

48 

78 

jellium 

4 

8 

12 

19 

25 

2 

full 

10 

29 

62 

114 

189 

jellium 

6 

15  a) 

24  b) 

39 

53 

3 

full 

17 

50 

111 

207 

347 

jellium 

9 

22 

38 

64 

90 

The  weighted  standard  deviations  of  the  fits  to  the  levels  in  methane  are 

a)  0.0056  cnf^  fit  to  levels  through  J=12. 

b)  0.0205  cm  ^  fit  to  levels  through  J=20. 


Table  6  Vibrational  mode  frequencies  for  C!!^  as  suggested  by  a  perturbation 
expansion  based  on  a  spherically  symmetric  unperturbed  state. 


Spherical  Symmetry  T,  Svmmetrv 

•  a 


a)  A  number  representation  is  used  in  which  the  number  of  vibrational  modes 
for  angular  momentum  £  is  written  as  n^  and  the  set  of  modes  by 

*£  "  <n0»nl»n2)* 

b)  In  the  overtone  modes  (i.e.  nj,  >  1) *  not  all  allowed  ^  values  occur 
because  of  the  identical  nature  of  the  particles  (phonons).  The  functions 
must  be  even  under  the  exchange  of  coordinates  of  identical  particles. 

c)  This  is  a  rough  estimate  of  the  mode  position  before  the  application  of  the 
Ta  perturbation.  The  indicated  accidental  degeneracy  between  the  various 

»s  not  symmetry  imposed  and  should  be  lifted  within  the  spherical 

approximation. 


d)  The  mode  designation  in  Td  symmetry  indicates,  the  number  of  Vi  modes  by  nv 
and  the  set  of  modes  by  it  =  [n  »  n\>  ’  nv>  *  n\>  ^  ’ 

vi  Vi  2  3  4 

e)  The  modes  positions  are  experimental  measurements  of  F2  symmetry  modes  as 
listed  in  Ref.  10  and  also  estimates  based  upon  the  multiples  of  the 
fundamental  frequencies  shown  in  brackets. 

f)  The  modes  involving  an  odd  number  of  £  =  1  vibrations  when  taken  in  Td 
symmetry  must  all  be  considered  as  a  direct  product  with  A«. 


